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ABSTRACT 

In  this  paper  we  investigate'  the  nonlinear  development  of  the  most  nnstal)le 
Gortler  mode  within  a  general  three-dimensional  l)onndary  layer  upon  a  suitably 
concave'  surface.  The'  structure  of  this  mode  was  first  identified  by  Denier,  Hall 
Se'ddeiugui  ( 1991)  whe)  de'monstrate'd  that  the  growth  rate  of  this  instability  is  0{G^  ) 
where  G  is  the'  Gehtle'r  numbe'r  (take'ii  to  be  large  lie're),  which  is  effectively  a  measure 
of  the'  curvature'  e)f  the  surfae-e.  Previe)tis  researclu's  have  described  the  fate  e)f  the' 
nie)st  unstable'  mode'  within  a  twe)-elimensional  bounelary  laye'r.  Dernier  Hall  (1992) 
elisceisse'el  the'  full}'  ne)nlinear  de'velopment  e)f  the  vente'x  in  this  case'  and  shejweel  that 
the'  neenline’arity  cause's  a  breakdown  eff  the  flenv  stnicture. 

The'  e'ffe-e-t  e)f  e-re)s.sfle)w  and  unste'adiness  eipeni  an  infinitesimal  unstable  meiele'  was 
elucielated  by  Bassenn  Hall  (1991).  They  dememstrated  that  crossfle)w  te'iids  te) 
stabilise'  the  me)st  unstable  Geirtle'r  mode,  and  for  certain  cre^ssflow/freque'ncy  e-ombi- 
nations  the'  Geirtler  mode  may  be'  maele  ne'utrally  stable.  These  vortex  cemfiguratiems 
naturally  leiul  the'inse’lve’s  te)  a  weakly  nonlinear  stability  analysis;  work  which  is  de- 
scribe'd  in  a  pre'vious  article  by  the  i)res<'nt  authors.  He're  we  e'xtend  the  ideas  e)f 
De-nie'r  and  Hall  (1992)  te)  the  thrc’e  dinK'nsie)nal  boundary  layer  pre)ble'm.  It  is  femnel 
tliat  the'  nume'iical  se)lutie)n  e)f  the  fully  ne)nhn('ar  eepiations  is  be'st  ce)nelucteel  using  a 
nie'the)el  which  is  e'sse-ntially  an  aelaption  e)f  that  utilised  ))y  De'iiie'r  anel  Hall  ( 1992).  The' 
infiue'iie-e'  e)f  e're)ssfie)W  and  unste'aeline'ss  upe)n  the  br<'akde)wn  e)f  the  fle)W  is  de'seribed. 
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§1  Introduction 


The  aim  of  this  article  is  to  further  our  understanding  of  the  effects  of  unsteadiness 
and  crossflow  upon  the  fully  nonlinear  development  of  unstable  Gdrtler  modes.  The 
initial  derivation  of  the  governing  equations  for  these  modes  was  given  l)y  Gdrtler 
(1940)  whose  results  were  modified  by  Hammerlin  (1956).  These  early  works  ignored 
nonparallel  efft'cts  present  due  to  boundary  layer  growth  and  Smith  (1955)  added 
terms  in  an  attempt  to  rectify  this  deficiency.  Until  recent  years  it  was  unclear  as  to 
tht'  importance  of  the  nonparalh'l  terms;  this  question  was  resolved  by  the  residts  of 
Hall  (1982  a,b,  1983).  In  this  series  of  papers  Hall  showed,  using  a  combination  of 
asymptotic  and  numerical  technitpies,  that  for  order  one  wavenumber  vortices  there  is 
no  unique  neutral  linear  stability  curve.  More  precisely,  the  stability  characteristics  of 
such  wavt'uumber  modes  are  entirely  dependent  upon  the  initial  form  and  location  of 
the  disturbance.  Howevt'r  for  small  wavelength  vortices  a  unique  ireutral  curve  does 
('xist  and  on  this  curv<'  the  vortex  wavenumber  k  is  O(G^)  where  G  is  the  (assumed 
large')  Gdrth'r  number.  For  hirther  reading  concerning  the  development  of  the  stability 
theory  for  Gdrtler  vortices  the  reader  is  referred  to  Hall  (1990). 

The  most  unstable  Gdrtler  mode  (i.e.  that  infinitesimally-sized  vortex  which  has 
th('  largest  growth  rate)  was  obtained  by  Denier,  Hall  Seddougui  (1991)  and  Tim- 
oshin  (1990).  By  considering  the  stability  properties  of  the  essentially  viscous  modes 
of  wavenumber  0(6'^  )  together  with  those  of  the  inviscid  modes  of  0(1)  wavelengths 
it  was  possibh'  to  identify  an  int<'rmediate  wavenumber  regime  in  which  the  vorte'x 
growth  rate  is  largest.  These  unstable  modes  are  confined  to  a  region  of  dejhh  of 
0(0“  s)  and  poss('ss  growth  rates  of  size  0(0®).  The  stability  properties  of  0(0^) 
wavenumber  vortices  are  deduced  by  solving  a  sixth-order  ordinary  differential  system 
with  appropriate  boundary  conditions  and  the  solution  of  this  system  reveals  that  the 
unicpie  most  unstable  mode  has  wavenumber  k  =  0.47CO5  and  growth  rate  0.3120^. 

To  dat('  there  has  be(’n  relatively  little  attention  paid  to  the  nonlinear  stability 
{)roperties  of  Gdrth'r  vortices.  Perhaps  the  first  work  was  performed  by  Aihara  (1976) 
who  attempted  to  describe  the  nonlinear  evolution  of  Gdrtler  vortices  using  parallel 
arguments.  Later  calculations  by  Hall  (1988)  remedied  these  defects  and  showed  that 
as  0(1 )  wavenumljer  modes  evolve  downstream  so  the  energy  of  the  flow  concentratt's 
its('lf  in  the  fundamental  and  mean  flow  correction.  This  suggests  far  downstream  the' 
flow  can  be  adequately  described  by  a  mean  field/first  harmonic  structure  and  such  a 
configuration  was  elucidated  for  short  wavelength  modes  using  both  weakly  nonlinear 
and  fidly  nonlinear  approaches  by  Hall  (1982  b)  and  Hall  Lakin  (1988)  respectively. 

.4  fuller  descrij)tion  of  the  nonlinear  stability  of  Gdrtler  modes  may  Ix'  found  in 
Denier  4;  Hall  (1992).  In  that  pai)er  the  authors  arg\ied  that  in  a  numlx'r  of  practical 
systc'ins,  espc'cially  where  significant  curvature  occurs  such  as  tlx'  case-  of  flow  ovt'r 
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turbine  blades,  one  would  expect  the  small  localized  surface  imperfections  may  well 
trigger  the  most  unstable  linear  Gortler  mode.  (This  conclusion  relies  on  the  result 
of  the  linear  receptivity  theory  given  by  Denier,  Hall  L  Seddoug\ii  (1991).)  This 
motivated  a  careful  study  of  nonlinear  evolution  of  the  most  unstable  mode  which  was 
tackled  in  the  following  way.  Denier  &  Hall  (1992)  took  an  arbitrary  form  of  initial 
disturbance  at  a  typical  streamwise  location,  say  .r  =  1.  By  integrating  the  linear 
form  of  the  vortex  equations  over  a  long  distance,  up  to  .r  =  101.  a  flow  profile  was 
obtained  which  is  dominated  by  the  most  unstable  mode  described  by  Denier  ct  al 
(1991).  From  x  —  101  onwards  the  full  nonlinear  equations  replaced  the  linear  ones 
and  the  most  unstable  vortex  mode  was  inarched  further.  Typically  8  or  16  harmonics 
of  the  fundamental  were  retained  during  this  calcidation.  It  was  found  that  at  a  critical 
jioint  the  flow  contains  a  region  of  reverse  flow  and  the  analysis  is  then  no  longer  valid. 
Denier  Hall  (1992)  interprete'd  this  breakdown  as  being  responsilile  for  the  vortices 
moving  away  from  the  wall  and  into  the  core  of  the  boundary  layer. 

The  effect  of  crossflow  and  unsteadiness  on  the  most  unstable  Gortler  mode  was 
discussed  by  Bassom  &  Hall  (1991).  The  primary  result  arising  from  this  work  was 
the  demonstration  that  a  relatively  small  crossflow  could  completely  stabilise  the  most 
unstable  mode.  Additionally,  by  allowing  for  vortex  unsteadiness,  it  was  shown  that 
suitable  combinations  of  crossflow,  vortex  frequency  and  wavenumber  could  lead  to 
neutrally  stable  coufigviratious.  .4  weakly  nonlinear  stability  analysis  pertaiiring  to 
such  configurations  was  conducted  by  Ba.ssom  4:  Otto  (1992)  who  derived  classical 
“Stuart -Watson'  (1960)  type  evolution  equations  for  near-neutral  modes.  They  con¬ 
cluded  that  the  weak-  nonlinearity  has  a  stabilising  effect  and  derived  equations  for 
the  supercritical  equilibrium  amplitudes. 

The  results  of  Denier  4:  Hall  (1992)  and  Bas.som  Hall  (1991)  provide  the  mo¬ 
tivation  for  the  current  study.  Within  a  two  dimensional  boundary  layer  the  effect  of 
nonlinearity  on  the  most  unstable  mode  tends  to  lead  to  a  finite-distance  breakdown 
wlu'reas  crcrssflow  aj)pears  to  stabili.se  the  flow.  With  these  two  mechanisms  tending 
to  have  ojipositf'  effects  it  is  clear  that  in  many  practical  situations,  where  three 
dimensionality  is  undoubtedly  imi)ortant,  it  is  of  great  interest  to  determine  which 
of  these  two  conflicting  behaviours  dominate.  We  attempt  to  answer  this  question  by 
considering  the  full  nonlinear  vortex  equations  and  employ  numerical  techniques  which 
an'  similar  to  those  u.sed  in  Denier  4:  Hall  (1992)  but  modified  in  certain  ways  (detailed 
lat('r).  These  improvements  significantly  speed  up  the  computations  and  allow  us  to 
olitain  a  gn'ater  range  of  results  than  tho.se  found  by  Denier  4:  Hall  (1992). 

The  structure  of  the  remaiiuh'r  of  this  article  is  as  follows:  in  section  2  the  fully 
nonlinear  cfiuations  are  dt’iived  and  a  brief  description  of  the  numerical  jnocedures 
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used  are  outlined  in  section  3.  Details  of  the  results  are  presented  in  section  4  and  in 
section  5  some  brief  conclusions  are  drawn. 


§2  Formulation 

As  in  Hall  (1985)  we  consider  a  boundary  layer  flowing  over  the  cylinder  y  =  0, 
—  oo  <  5  <  oo  where  the  5-axis  is  a  generator  of  the  cylinder,  y  measures  the  distance 
normal  to  the  surface  and  x  denotes  the  distance  along  the  surface.  The  Reynolds 
number  Re  and  Gortler  number  G  are  defined  by 

UqL  i 

/?.  =  —,  G  =  2R^S, 
u 

where  Ua  is  a  typical  flow  velocity  in  the  x-direction,  L  is  a  characteristic  streamwise 
lengthscale  and  v  is  the  kinematic  viscosity  of  the  fluid.  Furthermore  the  curvature  of 
the  cylinder  is  supposed  to  be  |xo  (x)  and  with  these  definitions  8  =  L/6,  where  h  is 
a  typical  radius  of  curvature  of  the  cylinder  (Gortler  vortices  are  typically  observed  in 
flows  over  concave  surfaces  which  corresponds  to  the  choice  Xo  >  0  ).  The  Reynolds 
number  is  supposed  to  be  large  whilst  8  is  sufficiently  small  so  that  as  ^  0  the 
parameter  G  is  fixed  and  is  of  order  one.  The  basic  three-dimensional  boundary  layer 
is  taken  to  be  of  the  form 

u  =  f7o  {u{X,Y),Rf^v{X,Y),R:h*w{X,Y))  (l  +  o(i?r^)) 

i 

where  X  =  x/L  and  Y'  =  yR?  jL  and  the  crossflow  parameter  A*  is  of  order  one. 

I 

It  is  convenient  to  define  the  scaled  .spanwise  coordinate  Z  =  zRJ  I L  and  let  T  be 
the  temporal  variable  scaled  on  L/Uq.  The  basic  velocity  profile  is  perturbed  by  the 
cpiantity 

(i7(T,x,y,z),R;"F(T,x,r,z),Rr%(r,x,F,z)) , 

_  i 

and  the  pressure  field  by  Re  P  (T,  X,Y,  Z).  Substituting  this  flow  forr.  into  the 
continuity  and  Navier-Stokes  equations  yields  the  .system 


Ux  +  VV  +  Wz  =  0, 


(2.1u) 


—Ut  +  Uyy  +  Uzz  —  uyV  —  uUx—uxU  —  ^^UY  -  \'wUz 

=  UUx  +  VUy  +  WUz, 

—  Vt  +  VVy  +  Yzz  —  GxuU  —  Py  —  TiVx  —  Y’xU  —  ^Vy  —  vyV  —  A*TeV'. 

=  UVx  +  FVV  +  +  ^\U\ 


(2.1/1) 


(2.1c) 
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-Wt  +  Wyy  +  Wzz  -Pz-  uWx  -  yWxU  -  Ftw  -  yVwY  - 

{2.1(1} 

=  UWx  +  VWy  +  WWz. 

—  1 

where  terms  of  relative  order  0{Rf.  ^ )  have  been  neglected.  It  is  worth  noting  at  this 
point  that  the  lineari.sed  system  studied  by  Bassom  Hall  (1991)  is  obtained  by  sc'tting 
the  right-hand  sides  of  (2.1  l)-d)  to  zero  whereas  the  nonlinear  equations  ('xaminerl 
by  Denier  &  Hall  (1992)  to  determine  the  development  of  steady  nonlinear  vortices  in 
two  dimensional  boundary  layers  can  be  retrieved  by  setting  A*  and  dr  ecpial  to  z('ro. 

We  now  invoke  the  scalings  jwoposed  by  Denier  e.i.  al.  (1991)  who  demonstrated 
that  in  high  Gdrtler  number  flows  the  most  unstable  vortices  have  0(G^  )  wavenumbers 
and  are  confined  to  a  layer  of  thickness  of  0{G~^)  adjacent  to  the  cylinder.  These 
modes  have  a  spatial  growth  rate  of  0{G^)  and  we  n.se  the  results  of  Bassom  Hall 
(1991)  who  illustrated  that  the  three  dimensionality  of  the  basic  flow  significantly 
affects  the  two  dimensional  stability  results  once  the  scaled  crossflow  parameter  A* 
becomes  0{G^).  Therefore  it  is  convenient  to  <lefine  the  0(1)  crossflow  parameter  A 
by 

A*  =  O^A.  (2.2a) 

To  reflect  the  fact  that  the  vortices  are  confined  to  a  re’gion  of  depth  0(G~^  )  adjacent 
to  the  cylinder  we  introduce  the  0(1)  boundary  layer  coordinat('  y  defined  by 

y  =  Gh\  (2.2/>) 


and  in  this  layer  the  basic  flow  may  be  expandefl  as  a  Maclaurin  series  of  the  form 


It  =  G  i//ii(A'),(/  +  2-0  +  j:G  S//.]3(A')/y'’  H - 

,r  =  G~^ft2\{X)y  +  ^G~^  fi22iX)}/  +  ^G~'^'  n2s{X)t/  +  •  •  ■ 


(2.2c) 

{2.2(1} 


To  determine  the  form  of  tin'  vortex  distvirbanc«*  we  appeal  to  the  findings  of  Bassom 
ik:  Otto  (1992)  who  identified  the  crucial  pertiirbation  size  at  which  th('  governing 
('(juations  b(>come  fully  nonlinear  (altlunigh  these  authors  made  no  attempt  to  solve 
these  fully  nonlinear  forms).  The  disturbance  forms  are  then 


=  c 

I  0  +  Ct  1  ]  T  G  t  2  "b  •  •  , 

r  =  G^  (Vn  +  G~U]  +  G-h2  +  ■•■), 

{2.3(1. h} 

W  -  1 

G*  ^  ^ n  0  -f  G  H  1  +  G  ^  H  2  4"  ■  •  ■ 

P  =  G^  (Po  +  G-^P, +G-^/^2+  -  ) 

(2.3c.  d) 

where  f  II,  In.  ITn.  Po.  U\ ....  etc.  are  all  functions  of  A’,  //.  th('  temporal  and  th<’  sjian- 
wise  variabh'.  It  is  now  convc'uient  to  implement  the  rc'sidt  of  Bassom  ic  Hall  (1991) 
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that  the  leading  order  behaviour  in  the  downstream  coordinate  is  purely  oscillatory 
and  we  do  this  by  introducing  the  coordinate  and  temporal  variations  given  by 


X  =  G-'^x,  Z  =  fSx  +  G-^z,  T  =  G-^t, 


(2.4a) 


where  (S  =  Ap2i//Mi-  This  then  implies  that  the  streamwise  and  spanwise  derivatives 
become  „  „  ^ 


The  desired  governing  equations  are  obtained  by  substituting  (2.2)-(2.4)  into  equations 
(2.1).  Leading  order  terms  in  the  momentum  equations  yield  that  Wo  =  I^Uo-  Next 
order  terms  in  equation  (2.1  b)  give  the  first  relationship 


a  d  d  d  \ 

V  +  a?  -  -  af  =  "“S.. 


(2.5o) 


where  the  precise  form  of  RHS|  is  given  presently  and  a  =  A//22  —  -A  second 

equation  is  derived  by  following  a  procedure  similar  to  that  described  by  Bassoni  &  Hall 
(1991).  By  considering  the  second  order  terms  in  the  y  and  5  momentum  equations 
(2.3),  eliminating  the  pressure  by  cross-differentiation  and  applying  the  continuity 
equation  it  is  a  routine  but  lengthy  task  to  obtain 

ydy^'^dz'^  2  dz  dT  dx )  \dy^  dz^ )  dz  ^ 

(2.56) 

where  again  we  shall  specify  RHS2  shortly.  The  forms  of  RHS2  includes  reference  to 
the  combination  H'l  —  /fCq  which  therefore  needs  to  be  expressed  in  terms  of  quantities 
with  subscript  zero.  This  is  best  accomplished  by  integrating  the  continuity  equation 
to  give 

T,.  .  f  \dUo  ,  dVo]  ,  , 


Wi  -  fiU]  =  4)0  (x 


dx  dy 


(2.5c) 


The  aim  of  this  article  is  to  consider  the  nonlinear  evolution  of  modes  which  are  periodic 
in  the  spanwise  direction,  with  fundamental  wavenumber  (\o/iii)  *  k  say,  and  periodic 
in  time.  Now  it  is  advantageous  to  introduce  scalings  first  proposed  by  Bassom  &  Hall 
(1991).  It  is  found  that  if  we  transform  according  to 

dy  (xohii)"  kdy,  dz  ->  (\o//ii) "  kd,,  dr  xh^u^x,  (2.6a  -  c) 


Uo  =  Xo"^^u^^^  ^o  =  Xo/<n^’  dr  ^  Xof^udr^ 
,  III  \  A  (/fii/t22  -  /<2I/'12) 

4>o  =  \of<u4), 


{2M-f) 


(2.6/7  -  6) 


o 


tlu'ii  the  transfonned  equations  (2.5  a,b)  are  independent  of  the  particular  boundary 
layer  under  consideration.  Therefore  the  ensuing  results  are  potentially  relevant  to 
a  wide  class  of  basic  configurations.  In  this  nondimensionalised  coordinate  system 
w('  restrict  ourselves  to  flow  quantities  the  fundamental  mode  having  spanwise  and 
tem])oral  variation  given  by  exjj  [?  (~  +  ^  re'presents  the  nondimensional 

fretpieiicy  of  the  fundamental  vortex  component. 

Upon  making  the  transformations  (2.6)  the  leading  order  vortex  ecpiations  (2.5 
a,]))  become 


2AdV  y  d'^V  I  I  1 

*  "  '  1  jl..)  £)_  ~  ^.3  Q~2  ~  f.2  Q^2  ~  0-0..  ~  ;.3  '  (--Ta) 


dt/ 


dzdy  dxdy 


k:^  k 

wh('r<’  th('  oj)erator  L  is  defined  according  to 


d^il'  \y^  d4'  1  dll'  y  dxl> 

^  “  BWf  ~  F 


(2.7A) 


(2.7c) 


and 


‘kV-lL 


dz 


d-U\ 


(2.7d) 


+  Vk  ^ 


^  ^  dy  J 


(2.7c) 


.01 


0. 


(£*•!)■ 


ou  .  ou 


{2.7f) 


To  com])l('t('  the  governing  ecpiations  it  is  neccvssary  to  determine  0.  This  is  achieved 
by  noticing  that  lUi  —  dt’i  satisfies 

OP 

I(ir,  -  dU, )  =  oiyUn  + 

Oz 

Idc'iitifying  the  comjionc'iit  of  this  ecpiation  independent  of  c  yields 

'2n 


0^  y  0 
“  F 


+  dz 


(2.7 g] 


G 


To  close  the  system  requires  specification  of  appropriate  boundary  conditions.  Clearly 
we  require  the  velocity  comjjonents  V  to  vanish  on  y  =  0  (and  by  continuity  so 
must  dV/dy).  Additionally  the  mean  flow  terms  =  0  on  y  =  0  and  in  ordc'r  that 
the  disturbance  b«'  confined  to  the  boundary  layer  we  flemand  that  the  streamwise 
velocity  U  tends  to  some  function  of  x  as  y  — »  co. 

We  notice  that  equations  (2.7)  are  the  appropriate  generalisations  of  those  solved 
by  Denier  and  Hall  for  the  most  unstable  noidinear  vortex  within  a  two  diiiK'nsional 
boundary  layer  (their  ecpiations  are  recovered  by  setting  A  =  ^  =  0  and  sc'tting  57  =0). 
This  allowed  us  to  compare  numerical  results  against  their  previously  published  ones 
as  a  check  of  our  numerical  methods. 

§3  Numerical  Methods 

The  methods  employed  to  solve  system  (2.7)  are  similar  in  essence  to  those'  use'd 
in  Denier  &  Hall  (1992).  However  it  was  found  necessary  to  introduce  a  numbe'r  of 
ame'iidments  to  their  code  in  order  to  speed  up  the  computations  which  in  turn  had 
the  benefit  of  allowing  an  increased  number  of  results  to  be  obtaineel. 

Denier  Hall  (1992)  solved  the  two-dimensional  counterparts  of  equations  (2.7) 
as  follows.  They  decomposed  each  of  the  flow  velocities  U,V,W  into  their  Fourier 
comjjonents  and  rewrote  the  governing  equations  in  terms  of  th('s('  components.  By 
utilising  a  scheme  based  upon  that  implemented  by  Hall  (1988)  they  obtained  finite- 
difference  ecpiations  for  the  mean  flow  and  harmonic  terms.  These  equations  contain 
only  one  streamwise  derivative  and  a  straightforwartl  method,  based  upon  solving  one 
tridiagonal  and  one  pentadiagonal  system,  may  be  used  to  march  the  solution  from 
one  streamwise  station  to  the  next.  For  more  details  of  the  practicalities  of  the  scheme 
the  reader  is  referred  to  Denier  Sc  Hall  (1992). 

One  marked  difference  betwet'ii  our  present  work  and  that  of  Denier  Sc  Hall  (1992) 
is  that  we  chose  to  calculate  the  nonlinear  terms  in  physical  space  rather  than  trans¬ 
form  space,  as  this  is  a  computationally  ‘cheaper’  policy  anrl  so  allowed  us  to  retain 
more  modes  in  our  calcidations.  This  necessitates  transforming  from  Fourit'r  spac(’ 
to  physical  space,  effecting  the  calcidations  and  extracting  the  Fourier  coefficients.  It 
was  decided  to  employ  Fast  Fourier  Transforms  to  do  this  which  has  th('  benefit  of 
reducing  the  cost  of  the  transformations  from  O(N^),  (the  cost  of  reducing  N  modes 
and  using  the  Fourier  transform  directly)  to  0(5iVlogiV  —  QN).  The  code  used  was 
based  on  the  original  Cooley  and  Tukey  (1965)  algorithm  and  was  thus  limited  to  .V 
lu'ing  an  integral  power  of  2. 

Two  other  changes  were  made  to  the  code  used  by  Denier  Hall  (1992).  As 
was  flescribed  in  tlu'  introduction,  these  authors  were  concerned  with  tlu’  devi'lojuiK’iit 
of  the  most  unstable  Gdrtler  mode  within  a  two  dimensional  boundary  layi’r.  To 
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achieve  this  most  unstable  mode  Denier  Sc  Hall  (1992)  took  an  arbitrary  disturljance 
and  then  marched  the  linear  stability  equations  for  a  long  distance  downstream  Viefore 
changing  to  the  nonlinear  system.  This  of  course  had  the  effect  of  ensuring  that  the 
most  unstable  component  of  the  original  disturl>ance  was  then  dominant  over  all  the 
otlu'r  components  so  that  when  the  fidl  nonlinear  equations  were  invoked  the  iiqjut 
vortex  flow  was  dominated  by  the  most  unstable  mode.  A  drawback  of  this  method 
is  that  the  lead  in  time  in  which  the  linear  mode  develops  proves  to  be  a  significant 
proportion  of  the  total  computation  time.  To  alleviate  this  difficidty  we  used  the 
program  described  by  Bassom  Sc  Otto  (1992)  to  compute  the  unstable  eigenfunction 
directly.  Tins  liad  the  effect  of  making  the  lead-in  time  to  the  nonlinear  equations 
redundant  and  tints  we  could  start  our  computations  of  the  full  nonlinear  system  almost 
immediately.  Substantial  reduction  in  the  computational  time  was  also  achieved  Ity 
using  a  stretched  grid  in  the  (/-direction  (normal  to  the  cylinder  surface)  as  opposed 
to  a  uniform  one,  and  the  scheme  eventually  chosen  for  this  jtrocess  was  taken  from 
Macaraeg,  Streett  Sc  Hussaini  (1988).  A  grid  is  required  which  encompasses  the  region 
Ix'tweeii  i/f,  =  0  and  ijp  =  where  the  subscript  p  denotes  the  physical  coordinate 

and  y,„„r  denot('s  some  outer  bound  at  which  the  a,symptotic  forms  of  the  solution  of 
the  system  (2.7)  are  supposed  to  have  been  attained.  A  notional  computational  grid 
0  <  <  1  is  introduced  which  is  related  to  pp  via 

ymnx  (^2  1)  l/r 

Up  —  TTjq  ^ 

(o  -  yt) 

with  the  vaUu'  (c2  —  1)  taking  values  between  0.2  and  1.  The  quantity  (c2  —  1)  rep- 
r(’sents  the  degree  of  stretching  between  the  large  and  small  steps  whilst  ej  controls 
the  rate  of  stretching,  1  <  C]  <  G.  For  most  of  the  calculations  reported  here  we  used 
ci  =  2.4  and  C’>  =  2  which  allows  100  grid  points  to  be  distributed  between  (/p  =  0 
and  t/p  —  50.  The  resulting  distribution  of  points  and  the  corresponding  step  lengths 
are  illustrated  in  figure  1.  Notice  that  as  y  y,nai  step  lengths  become  greater 
than  unity  which  in  a  finite  difference  scheme  may  induce  errors.  In  the  present  case 
gradients  of  the  dept'udeut  functions  in  this  regime  are  minuscule  and  so  this  possible 
difficulty  did  not  ari.se.  This  jiarticular  non  uniform  stretching  of  the  grid  points  al¬ 
lowed  a  twf’iity  fold  reduction  in  the  number  of  grid  points  over  the  regular  grid  used 
in  Bassom  Sc  Hall  ( 1991)  and  a  four  fold  improvement  over  the  piecewi.se  constant  step 
length  grid  used  in  Bassom  Sc  Otto  (1992). 

As  pn'viously  mentioned  the  solution  strategy  used  here  is  essentially  that  de¬ 
scribed  in  Hall  (1988).  Siij^pose  that  the  .solution  (Uo^Vo^L^Ot-  Oo)  is  known  at  some 
specifif'd  station.  ,r.  ami  suppose  further  that  a  guess  is  made  for  the  sohition  at  .r  -|-  f . 

noTdinear  terms  on  the  right  hand  sides  of  equations  (2.7  a,b) 
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were  calculated  using  these  guesses  and  the  system  solved  to  provide  updated  values 
These  updated  variations  were  then  put  into  the  nonlinear  terms 
and  the  process  repeated.  When  the  difference  between  successive  iterates  was  than 
some  norm  convergence  was  deemed  to  have  occurred  and  this  process  marched  onto 
the  next  step.  For  a  computational  grid  compo.sed  of  M  points  and  a  calcidation  with 
N  modes  retained  then  the  convergence  norm  usod  was 


A.-iV/2-l  j=M 

E  E 

k=-N/2  }=\ 


It'’; 


c'.l’l  +  iv;',;' 


(»+i) 


(»)  I 


where  denotes  the  iterate  of  the  A'**'  Fourier  component  of  the  flow  cpiantity 

U  evaluated  at  the  y-position. 

The  equations  (2.7  a,b)  may  be  discretised  for  any  particular  mode  into  the  forms 


+  2  T  1  T  frn  T  I’m  - 1  T  Cm  t^m~2  d"  n ,,,  —  Vmi  (n?  —  3,  •  ■  •  d/  1) 

(3,ln) 

and 

T  T  ’ni’^m  — 1  d"  Jm^'m  ~  Uni-  (/U  =  2,  •  •  •  A/  1)  (3.16^ 


This  system  was  solved  using  a  technique  outlined  in  Appendix  B  of  Bassom  tL’  Otto 
(1992).  This  solves  the  problem  (3.1)  in  one  sweep  rather  than  treating  the  pair  (3.1 
a,b)  as  distinct  pentadiagonal  and  tridiagonal  problems.  In  essence,  the  equations  are 
solved  by  performing  an  outward  forward  elimination  followed  by  a  back  substitution. 

Earlier  we  alluded  to  the  fact  that  significant  computational  saving  was  obtained 
by  utilising  Fast  Fourier  Transforms  in  preference  to  decomposing  the  flow  quantities 
into  their  respective  components.  Standard  procedures  were  implemented  so  that  to 
change  between  physical  and  transform  spaces  we  write 


i'j  = 


d  e  [i,iV], 


--N 


E 

7=1 


t’C 


k  e 


Y'  Y 


1]. 


where  Zj  =  27r(j  —  l)/iV  for  j  6  [1,jV]  find  where  quantities  with  tildes  are  in  the 
tran.sform  space  and  those  without  in  the  physical  space.  Denier  Kz  Hall  (1992)  were 
able  to  restrict  themselves  to  using  a  cosine  basis  <lue  to  the  nature  of  Gdrtler  vorticc's 
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with  two  (liiiu’iisional  bovuiclary  l<iy«'rs  hut  tlio  addition  of  ti'inporal  jx'iiodicity  and 
crossfiow  prevented  us  from  doing  *'lns. 

For  tli('  mtijority  of  tlu'  ealeulations  j>resented  Ih’Iow  the  nuinher  of  inodt's  retained 
was  1C  and  a  fairly  large'  step  use'd  ui  the  streainwise'  dirt'ction.  typieallv  f  =  0.01. 
De'iiit'i-  iV-  Hall  (1992)  found  in  tlu'ir  ealeulations  that  these'  paraine'tt'r  ehoiees  gave 
re'sults  to  within  graphical  aeeuraev.  Further  te-sting.  using  100  points  in  the  normal 
dire’ction  and  with  =  50  has  eonfirme'el  that  simihir  choice's  are'  satisfacteuy  feer 

the'  thre'e'  dime'nsional  computations  pe'rfeunu'el  he'ie'. 

^4  Results 

WV  Inive'  detade'd  the'  nume'iie'al  me'thoel  hv  whie'h  we'  m\'('stig<ite'd  the'  solution 
pro[)e'rtie's  of  syste'in  (2.7)  td*^hough  as  yt't  we  have  left  the  de  finition  of  the  amplitude' 
of  the  initial  vorte'X  unspe'C’fieel.  As  descrihe'd.  we'  initiated  our  comiuitatioiis  with  a 
multiple  of  the'  e'ige'ufunctions  of  the'  line'arise'el  v('r.sii)ns  of  system  (2.7).  For  a  spe'cified 
vorte'X  wave'r.eunhe’r  k\  fre'epu'iie-y  Q.  anel  cre>ssflow  A  the  method  outlined  by  Bassom 
A’  Otto  (1992)  was  use-el  to  e-omi)ute  the  e'orre'sponeling  linearise'd  growth  rates  and 
the'  re'spe'e'tive'  e'ige'ufunctions  normalise'el  so  that  the  ene'rgy  de'fine'd  In¬ 


is  e-eiual  to  A'.  We'  refe'r  to  A  as  the'  amj)lituel<'  of  the*  initial  e'ondition. 

Our  first  e-alcul<itions  we’re'  pe’iforme’el  primarily  as  a  verifie’ation  of  our  e-ode  against 
the’  e’st<d)lishe’d  le'sults  of  De’iiie'r  A;  Hall  (1992).  fhe'nce'  we’  e-onsiele’red  the  case  of  a 
jeure’ly  t wo-dime’nsie)nal  be)unel<-iry  laye’r  and  ste’aely  ve)rtice's  ( A  =  If  =  0).  De’iiier  A: 
Hall  (1992)  founel  tlnit  as  the'  \-e)rte'X  e'vedve's  ne)nline’arly  elownstre’am  the’n  at  some' 
point  the-ir  coinpeit.-itienis  broke’  elown  it.  a  singedarity.  Care'ful  inve'stigatieen  re’ve’ale'el 
that  at  this  location  the'  skin  frictie)n  was  e)f  the'  fieew  vanislu's  anel  the'ii  ;uiv  sohition 
se  he’ine’  whie-h  re’lie’s  upon  a  mare-hing  te'e-hideiue’  be'e-ome's  invalidate’el.  The’  vanishing  eif 
the’  skin  frictiein  was  inte’rpre’te’d  as  be’ing  iiulicative'  eif  the’  veirtie-e’s  breaking  awav  freim 
the'  wall  tuiel  moving  intei  the’  e-ore'  eif  the'  bounelary  Itiye'r.  Our  re’sults  for  varieius  initial 
;ini])lituel;’s  A  at  ’  inelie-ate-el  eiii  figure’  2  whe'ie'  we'  sheiw  the’  leie'atieui  eif  the'  bre’akele^wn 
jieeint  .I'h  ;is  a  funi'tion  eif  A  feir  the*  me)st  unstable’  Gdrtle’C  moele’  with  veirte'x  wave'nundie'r 
k  =  0.470  (the-  disi  :bane-e'  was  intrexhice'el  at  .r  =  1  sei  the'  distane'e'  trave’lle’el  liv  the 
pe'iturbation  be’feue  bre-akelenvii  is  .I’h  —  1).  Not  surprisingly,  the  bre’akelenvn  leM-atieui 
.(■(,  is  a  meeneeteuie'  ele’e-re'asing  functieni  oi  the  amplitude’  A  anel  as  A  — e  0  se)  Xh  -+  oc 
euiel  the’  line'Jir  probli’in  is  re’trie’ve'fl.  Fe^r  all  the'  eithi'r  e-alcuhitieins  re'peirte’el  upeui  he’re’ 
a  simil.-ir  tre’iid  is  eib  e-rve-el  so  that  in  all  case’s  we-  e’hexsc  A  =  0.2:  this  si'le'ctieui  was 
made’  pure’ly  for  ihustr<-iti\e'  purpeise's  anel  has  nei  sjie'cied  signihe-ane'e'  whatse)e'^-e'r. 
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Ill  figure"  3  w('  demonstrate  another  feature  of  the  two  elimensional  results.  Figure 
3a  indi.'ates  tlie  depende'iiee  of  the  linear  vortex  growth  rate*  d,.  as  a  function  of  the’ 
vorte'x  wavenumher  k  for  the'  first  two  inodes.  As  is  now  well  documented  by  Denier 
rt  al.  (1991),  De’uier  Hall  (1992)  the  most  unstable  mode  has  growth  rate'  = 
0.312  at  k  =  0.476  and  lir  — ^  0  as  A:  — >  0  eir  k  oo.  We  note  however,  that  feir  a 
significant  range  of  scaled  wavenumbers,  roughly  0.3  <  k  <  0.8,  the  vortvx  growth 
rate'  is  not  much  reduced  from  the  maximum  in  as  much  as  it  is  at  least  90%  eif 
the  maximum.  This  will  have  important  consequences  for  later  results.  Figure  3b 
ilhist rate's  breakdenvn  point  .r^  as  a  function  of  k  for  the  first  two  modes.  As  is  to  be 
e'xpe'e'te'd  for  fixed  initial  ainjilitude  A  and  fixed  wavenumlier  the  first  mode  always 
bre'aks  away  from  the  wall  before  the  second  one  and  this  trend  continues  ftir  higher 
mofle's.  However,  it  is  also  observed  that  for  the  primary  mode,  Xb  is  a  monotone' 
de'e-re'asing  function  of  k.  Therefore  for  any  fixed  initial  amplitude  it  is  not  the  most 
unstalile  moele  which  breaks  up  first.  Indeed  we  previously  remarked  that  for  a  whole' 
range'  of  wavenumbers  surrounding  the  most  unstable  value  linear  growth  rates  are  not 
te)e)  diffe'ie'iit  from  the  largest  growth  rate.  This  has  potentially  imiiortant  consequene'e's 
fe)r  a  numljer  of  practical  flows  as  it  demonstrates  that  the  breakdown  properties  of  the 
flow  are'  crucially  de'pe'ndent  on  the  nature  of  the  physical  characteristics  of  the  flow 
anel  are'  se'iisitive'  to  the  nature  of  the'  evolutiem  of  the  flow.  More  i)recisely,  suppose 
that  the  vortex  motion  starts  with  e'xtremely  small  amplit\ide.  Then  one  would  expect 
there  to  be'  a  large  distance  over  which  the  motion  ilevelops  essentially  in  a  linear 
manner.  During  this  time'  the  me)st  unstable  linear  mode  would  overwhelm  modes  of 
othe'r  wavele'ijgths  and  one-e  the'  vortice's  had  grown  sufficiently  so  that  nonlinearity  is 
important  the  fle)W  behaviour  wotdd  be  dominated  by  that  of  the  most  unstable  mode. 
On  the'  other  hand  if  the  initial  vortex  amplitude  is  not  tiny,  nonlinear  effect  are  like-ly 
te)  be'  important  a  relatively  short  distance  downstream  by  which  stage  it  is  unlikely 
that  the  most  unstable'  mode  dominaM's  the  others.  In  this  case  the  behaviour  of  the 
most  unstable  vortex  is  not  likely  te)  dictate  the  properties  of  the  breakdown  of  the 
fle)w. 

We  turn  now  to  consider  cases  with  no:,  zero  crossflows.  As  noted  by  Bassom  <1,: 
Hall  (1991)  we  may  restrict  our  attention  to  case’s  in  which  the  crossflow  parameter 
A  >  0  siiice  by  suitably  transforming  the  system  (2.7)  we  can  relate  flows  with  A  <  0 
to  appre)priate'  counterparts  with  A  >  0.  In  figure  4  we  recall  the  results  of  increasing 
crossflow  on  the  linear,  stationary  vortex  mode.  As  discussed  in  detail  by  Bassom  & 
Hall  ( 1991)  the  e’ffe’ct  cT  crejssflow  on  linear  vortex  structures  of  wavelenigth  0(G“  s  )  is 
primarily  a  stabilising  one.  Indeed  figure  4  illustrates  that  when  A  >  0.410  the  vort('x 
niorle  is  stal)ilis('d  for  all  wavenumbers  in  the  0(G^ )  regime. 
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Figures  5  ai:  b  show  the  variations  of  the  inaxiniuin  growth  rate  finfl 

the  rorresponding  vortex  wavenunil)er  k,nar  for  increasing  crossflow  paranieti'r  A.  As 
already  noted,  increasing  A  lowers  and  it  is  the  inode  with  wavenumber  k  ~ 

0.407  which  is  the  last  to  be  stabilised.  Further  is  not  a  monotone  function  of  A 
as  might  have  been  anticipated.  Figure  5b,  as  well  as  .showing  indicates  the 

growth  rate  of  the  mode  with  wavenumber  k  =  0.476  (the  value  of  k„tar  for  the  case  of 
zero  crossflow).  It  can  be  seen  that  the  difference  in  growth  rates  of  the  most  unstable 
mode  and  that  with  k  =  0.476  is  very  small  over  the  whole  range  of  A  considered. 
This  reinforces  the  earlier  comment  that  there  is  a  significant  wavenumber  regime  over 
which  the  linear  growth  rate  of  the  vortex  mode  is  almost  constant. 

For  each  crossflow  A  we  integrated  the  appropriate  most  unstable  linear  eigenmode 
with  initial  amplitude  A  =  0.2  until  breakdown  occurred,  and  then  repeated  the 
experiment  with  vortex  wavenumber  k  =  0.476.  The  results  are  summarised  in  figure 
5c.  It  is  seen  that  as  well  as  exercising  a  stabilising  influence  on  the  linear  mode, 
increasing  the  crossflow  tends  to  have  a  similar  effect  on  the  nonlinear  modes.  Clearly 
th('  greater  A,  the  greater  the  delay  downstream  before  the  initially  most  unstable 
mode  breaks  away  from  the  wall.  Additionally,  for  the  same  initial  amplitude  A  and 
crossflow  A  the  mode  with  k  =  0.476  breaks  uj)  before  the  linearly  most  unstable  mode. 
Further  investigations  have  suggested  that  for  all  choices  of  A  and  vortex  frequency 
Q  then  if  two  linear  modes  of  different  modes  kj  >  k2  are  such  that  if  initially  have 
the  sanu'  amplitude  then  if  they  are  marched  downstream  then  it  is  the  one  with  the 
higher  wavenumber  which  breaks  down  first. 

Finally,  we  discuss  in  more  detail  the  influence  of  unsteadiness  on  our  findings.  In 
their  study  Bassom  &  Hall  (1991)  made  some  comments  concerning  the  properties  of 
time  dependent  linearised  vortices.  For  the  majority  of  their  work  these  authors  were 
jirimarily  interested  in  examining  neutrally  stable  modes  although  they  did  compute 
a  few  non- neutral  ones  (see  their  figure  16).  Bas.som  &  Otto  (1992)  showed  that  for 
fixed  wavenumber  k  then  as  the  frequency  Q  of  the  mode  increases  so  the  crossflow 
reipiired  to  maintain  neutral  stability  grows.  In  particular  it  was  shown  that  for  a 
non  dimensional  vortex  frequency  the  stability  properties  of  the  vortex  are  sensitive 
to  the  sign  of  Cl.  In  figures  6  a&  b  we  illustrate  a  facet  of  this  sensitivity.  For  each 
frequency  fl  and  crossflow  parameter  A  we  show  the  wavenumber  of  the  most  unstable 
linear  mode,  kmai,  together  with  the  growth  rate  of  that  mode  <  0,  it 

it  is  observed  that  as  A  increases  from  zero  so  k^ar  decreases  whereas  for  >  0  this 
is  not  the  case.  Correspondingly,  when  Cl  <  d  the  growth  rate  of  the  most  unstable 
morh'  is  monotone  decreasing  in  A  whereas  when  >  0  then  as  A  increases  from  zero  so 
tliere  is  a  crossflow  range  over  which  increases.  This  increase  is  not  indefinite 

however  and  there  is  a  critical  A.  dependent  upon  Cl,  after  which  the  growth  rate 
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decreases.  We  also  note  that  if  Qi  >  ^^2  >  0  then  max^  /^r  <  niax^  ^$72; 


so  that  over  all  Cl  and  all  A  >  0  the  mode  with  the  greatest  growth  rate  is  stationary 
and  exists  in  a  two-dimensional  boundary  layer. 

The  breakdown  characteristics  of  unsteady  flows  are  described  by  figure  6c.  For 
each  frequency  Cl  and  A  we  marched  the  linearly  most  unstable  mode  of  amplitude 
A  =  0.2  from  .t  =  1.  We  can  observe  the  somewhat  conflicting  roles  played  by 


crossflow  and  frequency. 

In  the  main,  for  a  prescribed  Cl  increasing  A  delay.s  breakdown  whereas  for  fixed 
A  and  increasing  Cl  this  phenomenon  is  enhanced.  Notice,  however,  one  important 
feature  which  runs  against  this  general  rule.  For  positive  Cl  then  a  small  to  moderate 
crossflow  actually  tends  to  promote  breakdown  although  larger  crossflows  do  reverse 
this  effect.  An  attempt  was  made  to  verify  this  trend  by  considering  larger  values  of 
Cl  than  those  illustrated  in  figure  6.  However  problems  were  encountered  as  Cl  grew 
and  these  difficulties  can  be  attributed  to  a  number  of  causes.  Following  on  from  work 
elucidated  in  Bassom  &  Hall  (1991)  it  is  the  case  that  for  small  crossflow  the  most 
unstable  linear  mode  first  has  a  small  wavenumber  relative  to  the  implied  scaling.  As  A 
grows  then  the  most  unstable  mode  corresponds  to  an  eigenfunction  that  moves  away 
from  the  wall  at  y  =  0.  At  Cl  greater  than  about  2  this  movement  occurs  quickly  for 
small  changes  in  A  so  that  for  quite  moderate  values  of  A  the  eiegnmode  is  far  removed 
from  the  boundary.  As  found  both  by  Bassom  &:  Hall  (1991)  and  Bassom  Sc  Otto 
(1992)  the  numerical  solution  of  the  governing  equations  becomes  noir-trivial  as  the 
vortex  moves  out  since  boundary  conditions  need  to  be  imposed  at  the  wall.  Clearly 
for  modes  concentrated  away  from  the  wall  large  changes  in  /i,  A  or  Cl  can  lead  to 
almost  imperceptible  changes  in  the  values  of  the  eigenfunctions  at  the  wall  and  thus 
reliable  numerical  convergence  is  rendered  very  difficult.  However,  our  limited  further 
computations  for  >  1.5  are  in  agreement  with  the  general  behaviour  described  above. 


§5  Conclusions  &:  Discussion 

In  this  work  we  have  detailed  the  nonlinear  spatial  evolution  of  unstable  Gortler 
modes  in  a  three-dimensional  boundary  layer.  In  particular,  the  roles  played  by  vortex 
wavenumber,  frequency,  and  the  crossflow  component  of  the  underlying  base  flow  have 
been  described.  We  feel  that  of  particular  importance  is  our  finding  that  (  all  other 
factors  being  equal)  of  two  modes  of  wavenumber  within  the  O(G^)  regime  the  one 
v/ith  the  smaller  wavelength  will  be  the  first  to  breakdown.  This  then  suggests  that  in 
practical  situations  it  may  not  be  the  most  unstable  linear  mode  which  is  of  ultimate 
importance. 

In  many  cases  the  relevant  calculation  to  describe  the  lireakdown  of  a  flow  is  oiu' 
of  a  receptivity  type.  In  this  scenario  small  disturbances  within  the  boundary  layer 
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or  oil  the  wall  of  the  cylinder  can  trigger  Gortler  vortices  and  the  pri'cise  method 
of  this  triggering  frequently  excites  modes  of  a  jireferred  wavenumlier.  If  this  occurs 
in  practice  then  our  calculations  provide  a  descrijition  of  the  evolution  of  the  mod(\ 
Conversely,  if  a  range  of  wavenumber  modes  is  excited  two  eventualities  would  seem  to 
be  possible.  First,  suppose  that  the  initial  disturbance  is  very  small.  Then  it  is  to  lie 
expected  that  the  perturliation  travels  a  long  way  downstream  before  nonlinearity  has 
significant  effect,  the  most  unstalile  linear  mode  will  be  dominant  before  this  point  and 
its  evolution  characteristics  will  essentially  describe  that  of  the  whole  flow.  Second, 
suppose  that  the  initial  perturbation  is  not  so  small.  Our  results  summarised  by  figure 
3  have  shown  that  although  there  is  a  unique  most  unstable  mode  for  each  crossflow  A 
and  frequency  Q  vortices  with  wavenumbers  in  a  fairly-large  region  surrounding  that 
of  the  most  unstable  mode  have  growth  rates  not  very  different  iTom  the  maximum. 
Therefore,  by  the  time  nonlinearity  is  significant  it  is  not  clear  that  the  most  unstable 
mode  would  be  dominant  and  the  breakdown  characteristics  of  the  whole  flow  woidd 
involve  calculations  more  involved  than  those  reported  here.  However  we  have  shown 
that  for  given  mode  amplitude  it  is  the  higher-wavenumber  modes  which  appear  to 
))reakdown  first  so  that  these  components  of  a  spectrum  of  excited  modes  may  well 
prove  to  be  the  important  ones. 

Denier  Hall  (1992)  showed  that  when  their  calculations  for  nonlinear  modes  in 
two-  dimensional  boundary  layers  terminated,  this  corresponds  to  the  skin  friction  of 
the  flow  vanishing  at  some  point.  Once  this  happens  marching  schemes  as  used  both 
here  and  in  Denier  Sc  Hall  (1992)  cannot  be  continued.  We  confirm  this  finding  for 
our  three  dimensional  cases  as  well  but  we  also  observed  that  before  the  skin  friction 
vanishes  the  velocity  profile  develop  inflection  points  at  positions  away  from  the  wall. 
The  a])pearance  of  those  inflection  points  suggest  that  the  flow  will  become  susceptible 
to  Rayleigh  vvaves  which  would  give  an  alternative  route  to  the  idtiniate  breakdown. 
The  analysis  of  these  modes  would  be  of  interest. 

Finally,  we  recall  that  all  our  calculations  have  been  concerned  with  considering 
the  evolution  of  perturbations  of  a  siiecifled  wavenumber.  Of  course  in  some  situations 
a  spectrum  of  modes  may  well  be  present.  We  have  identified  situations  in  which  we 
might  expect  one  mode  to  dominate  the  others  before  nonlinearity  sets  in  but  in  the 
other  cases  calcidations  would  be  needed  which  account  for  an  initial  perturbation 
which  contains  a  number  of  modes.  The  development  of  a  code  to  perform  such  calcu¬ 
lations  might  well  be  formidable  but  it  would  give  the  definitive  theoretical  description 
of  nonlinear  Gortler  vortex  behaviour  in  three-dimensional  boundary  layers. 
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a)  Linear  vortex  growth  rates  iir  as  a  fiiiirtion  of  wavenumher  k  for  the  two  most 
flangerons  modes  in  a  two  rlimensional  ])ouiKlary  layer. 

h)  Breakdown  point  .rj,  as  a  fnnction  of  wavemimher  k  for  the  modes  in  figure  3a  with 
assnmefl  initial  amplitude  A  =  0.2.  The  symbols  denote  rorresponding  points  on  the 
curves  of  a)  and  b). 
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Figure  4 


A  increasing 


Figure  4;  Growth  rate  of  /ir  of  the  most  dangerous  stationary  vortex  modes  within 
a  boundary  layer  with  increasing  crossflow  A.  The  curves  correspond  to  A  varying 
between  0.005  and  0.410  in  9  equal  steps.  The  indicated  points  show  the  locus  of  the 
wavenumber  of  the  most  unstable  mode  as  the  crossflow  varies. 


a)  Wavriumilx’r  of  the  most  unstable  nonstationary  modes  and  1))  corrf'spoiirline; 
growth  rat('s  3r  as  fnnetions  of  A.  Here  we  have  eonsidered  Cl  =  -1, 0.5.  0.  0.5.  1.  1.5. 
e)  Breakdown  point  .rj,  as  a  fnnrtion  of  crossflow  A  for  the  modes  with  frecpu'iicic's  Cl 
and  wavennnrbers  k  as  given  in  figure  6a)  and  initial  amplitude  A  =  0.2. 
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